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ABSTRACT 

We present Keck laser guide star adaptive optics imaging and aperture masking observations of 
the M8+L7 binary LHS 2397aAB. Together with archival HST, Gemini-North, and VLT data, our 
observations span 11.8 years of the binary's 14.2-year orbital period. We determine a total dynamical 



mass of 0.146_q'q]^3 Mq (153_]^4 Mjup). Using the combined observational constraints of the total 
mass and individual luminosities, the Tucson (Lyon) evolutionary models give an age for the system 
of l.Sl^Q g Gyr (1.8lQ'g Gyr), which is consistent with its space motion based on a comparison to 



the Besangon Galactic structure model. We also use these models to determine the mass ratio, giv- 



ing individual masses of 0.0839 
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Q Mq) for LHS 2397aB. Because LHS 2397aB is very close to the theoretical mass-limit of 
lithium burning, which remains untested by dynamical masses, measuring its lithium depletion would 
uniquely test substcUar models. We estimate a spectral type of L7±l for LHS 2397aB, making it the 
first L/T transition object with a dynamical mass determination. This enables a precise estimate of 
its effective temperature from Tucson (Lyon) evolutionary models of 1450±40 K (1430±40 K), which 
is 200 K higher than estimates for young late-L companions but consistent with older late-L field 
dwarfs, supporting the idea that the temperature of the L/T transition is surface gravity dependent. 
Comparing our temperature estimate for LHS 2397aB to those derived from spectral synthesis mod- 
eling for similar objects reveals consistency between evolutionary and atmospheric models at the L/T 
transition, despite the currently limited understanding of this phase of substellar evolution. Future 
dynamical masses for L/T binaries spanning a range of surface gravity, age, and mass will provide the 
next critical tests of substellar models at the L/T transition. 

Subject headings: binaries: general, close — stars: brown dwarfs — infrared: stars — techniques: high 
angular resolution 



1. INTRODUCTION 

In the field population, the warmest brown dwarfs and 
the lowest mass stars are in many ways quite different. 
Beyond the fact that brown dwarfs by definition do not 

* Some of the data presented herein were obtained at the W.M. 
Keck Observatory, which is operated as a scientific partnership 
among the California Institute of Technology, the University of 
California, and the National Aeronautics and Space Administra- 
tion. The Observatory was made possible by the generous financial 
support of the W.M. Keck Foundation. 

t Based partly on observations made with the NASA/ESA Hub- 
ble Space Telescope, obtained from the data archive at the Space 
Telescope Institute. STScI is operated by the association of Uni- 
versities for Research in Astronomy, Inc. under the NASA contract 
NAS 5-26555. 

Based partly on observations obtained under program ID GN- 
2002A-Q-6 at the Gemini Observatory, which is operated by the 
Association of Universities for Research in Astronomy, Inc., un- 
der a cooperative agreement with the NSF on behalf of the Gem- 
ini partnership: the National Science Foundation (United States), 
the Science and Technology Facilities Council (United Kingdom), 
the National Research Council (Canada), CONIC YT (Chile), the 
Australian Research Council (Australia), Ministerio da Ciencia e 
Tecnologia (Brazil) and SECYT (Argentina). 

^ Based partly on observations made with ESO Telescopes at the 
Paranal Observatory under program IDs 071.C-0716 and 076. C- 
0139. 

^ Alfred P. Sloan Research Fellow 



sustain internal energy generation through hydrogen fu- 
sion, while very low-mass stars do, the effective temper- 
atures of these two classes of objects span a range of 
wlSOO K and more than an order of magnitude in lu- 
minosity. However, they share one important attribute, 
which is the presence of dust clouds in their atmospheres. 
Below effective temperatures of ^2800 K, photospheres 
are cool enough for particulate matter to form (e.g., Tsuji 
et al. 1996; Jones & Tsuji 1997), and this dust persists for 
objects as cool as wl300 K (e.g., Chabrier et al. 2000), 
at which temperature the dust begins to settle below the 
photosphere. Thus, the formation and sedimentation of 
dust clouds are among the key processes driving the at- 
mospheric physics of objects in this temperature range, 
which also includes the warmest extrasolar planets such 
as transiting hot Jupiters and young giant planets (e.g.. 
Barman 2008; Marois et al. 2008). 

In addition to providing excellent laboratories for ul- 
tracool atmospheric physics, very low- mass stars and 
brown dwarfs are also intrinsically interesting as the low- 
est mass products of star formation. In particular, the 
properties of binary systems spanning a wide range of 
stellar masses can provide discriminating tests of various 
star formation models. While the binary properties of 
low-mass stars (FGKM spectral types) have been exten- 
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sively studied (e.g., Duquennoy & Mayor 1991; Fischer 
& Marcy 1992), such studies have only become possible 
within the past decade for brown dwarfs and stars at the 
bottom of the main sequence. In fact, such studies are 
generally hampered by the fact that brown dwarfs follow 
a mass-luminosity-age relation, rather than the simpler 
mass luminosity relation for main-sequence stars. For 
example, a latc-M field dwarf may be a young brown 
dwarf or a star at the bottom of the main sequence. Bi- 
naries with well determined orbits arc extremely valu- 
able as they break this mass age degeneracy by provid- 
ing dynamical mass estimates, though very few mass de- 
terminations are available to date (e.g., Zapatero Osorio 
et al. 2004; Dupuy et al. 2009). In fact, mass measure- 
ments from binaries typically offer much stronger con- 
straints on substellar models than other commonly mea- 
sured parameters (e.g., age and compositon) given that 
their higher precision, hence the concept of using bina- 
ries as "mass benchmarks" (Liu et al. 2008). In addi- 
tion, binary orbits provide a wealth of other information 
such as the true (not projected) semimajor axis and the 
eccentricity, both of which have different predicted dis- 
tributions under different formation scenarios (e.g., Bate 
2009; Fisher 2004). 

One particularly striking result of the many surveys 
that have searched for brown dwarf binaries is the 
paucity of intermediate mass ratio binaries. Although 
such surveys have typically been sensitive to mass ra- 
tios as small as Mb/Ma « 0.6, the mass ratio distribu- 
tion drops sharply below unity, with > 50% of binaries 
having estimated mass ratios of > 0.9 (Burgasser et al. 
2007). This is in stark contrast to solar- type binaries 
whose mass ratio distribution peaks at «0.4 (Duquennoy 
& Mayor 1991), suggesting a different formation process 
for these two classes of objects. Although very unequal 
mass (or spectral type) brown dwarf binaries are rare, 
they provide the strongest coevality tests of substellar 
theoretical models because they offer the most leverage 
on constraining model isochrones. 

Freed et al. (2003) discovered a faint companion to 
LHS 2397a^ using the curvature AO system Hokupa'a 
at Gemini-North, revealing one of the most unequal- 
flux ultracool binaries known (Aif = 2.77±0.10 mag). 
LHS 2397a (M8) is an Ha flare star with a distance mea- 
sured by Monet et al. (1992) to be 14.3±0.4 pc. Freed 
et al. (2003) estimated the spectral type of LHS 2397aB 
from near-infrared photometry to be L7.5±l. Thus, 
LHS 2397aAB represents a rare pairing of two objects at 
the extreme ends of the effective temperature range over 
which dust clouds play an important role in atmospheric 
physics. LHS 2397aA is just cool enough that dust is 
likely to persist in its atmosphere, while LHS 2397aB 
has nearly cooled to the point of becoming a dust-free 
T dwarf. 

Based on Keck laser guide star adaptive optics 
(LGS AO) imaging and aperture masking observations 
from our ongoing orbital monitoring program of ultra- 
cool binaries, we present here a dynamical mass for 
LHS 2397aAB. Combining our Keck data with archival 
Hubble Space Telescope (HST), Very Large Telescope 
(VLT), and Gemini- North Telescope images, we mea- 

2 LHS 2397a is the star listed between LHS 2397 and LHS 2398 
in the LHS Catalog (Luyten 1979). 



sure a total mass of 0.146±0.014 M0, with the dominant 
source of uncertainty being the 3.0% error in the par- 
allax (which translates into a 9.0% error in the mass). 
This is the first mass measurement for a binary contain- 
ing an L/T transition brown dwarf, and we use the direct 
mass measurement along with the measured luminosities 
of both components to test evolutionary and atmospheric 
models in this temperature regime (f»1400 K) for the first 
time. 

2. OBSERVATIONS 
2.1. HST/WFPC2-PC1 

We retrieved HST archival images of LHS 2397aAB 
obtained with the WFPC2 Planetary Camera (PCI) on 
UT 1997 April 12 (GO-6345, PI Kirkpatrick). The data 
comprise four F814W images with exposure times of 2 s, 
1000 s, 1000 s, and 300 s, taken in that order in a single 
orbit. The companion LHS 2397aB is only seen in the 
latter three images, and the primary star LHS 2397aA is 
saturated in all but the first image. We used TinyTim 
(Krist 1995) to generate point-spread-function (PSF) 
models which were fit to the data in a similar fashion to 
our previous work (Liu et al. 2008; Dupuy et al. 2009), 
thus allowing us to determine the positions and fluxes 
of LHS 2397aAB in each image. Saturated pixels were 
masked in our PSF-fitting routine. The F814W flux ratio 
of the binary (4.18±0.08 mag) was determined by com- 
paring the mean peak flux (i.e., the model PSF's nor- 
malization constant) of LHS 2397aB in the latter three 
images to the peak flux of LHS 2397aA in the flrst image. 

Since positional information is carried almost entirely 
in the core of the PSF, especially for these undersampled 
WFPC2 images, we could not determine precise posi- 
tions for LHS 2397aA in the three images in which its 
core was saturated. Therefore, in order to measure the 
binary separaton and P. A., we had to somehow tie the 
latter three measurements of LHS 2397aB to the first 
shallow image of LHS 2397aA. This was accomplished 
by fitting a first-order polynomial to the positions of 
LHS 2397aB, in order to account for any telescope drift 
that occurred during the observations, and extrapolat- 
ing this fit to the unsaturated image of LHS 2397aA. 
We found a drift in the (x, y) position of LHS 2397aB of 
(0.82±0.03, 0.26±0.01) pix/hr over the 0.7 hr duration 
of the observations. After correcting for this drift and 
using a pixel scale of 45.54±0.01 mas/pix,^ we found a 
separation of 274±4 mas, and a position angle (P.A.) of 
87.3±0.8° (Table 2). 

To determine the imcertainties of our position and flux 
measurements we simulatcid the data using images of 
single stars, allowing us to assess systematic errors in 
our fitting procedure. We built a library of single stars 
from other programs that targeted brown dwarfs with 
ffST/WFPC2-PCl (GO-8563, PI Kirkpatrick; GO-8581, 
PI Reid; and GO-8146, PI Reid). We only used stars 
that had equivalent or higher signal-to-noise (S/N) than 
the science data so that we could degrade the S/N of 
the library images to match the science data. We also 
restricted ourselves to observations consisting of two or 
more images to allow robust rejection of cosmic rays. 

^ This is the quoted pixel scale from the WFPC2 Instrument 
Handbook for Cycle 13, which is consistent with other values in 
the literature (e.g., see the discussion by Liu et al. 2008). 
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Fig. 1. — HST, Gemini, VLT and Keck images of LHS 2397aAB shown chronologically by column. We do not rotate 
the HST data so that north is up in order to preserve the somewhat undersampled nature of the WFPC2 data. The 
Airy ring of the Keck PSF is visible in some Keck images. All images are shown on the same scale, I'/O on a side, using 
a square-root stretch for the grayscale images. Pixels within a 75 mas radius of the primary have been divided by 10 

to more clearly show the companion. Contours arc drawn at 0.62, 0.38, 0.24, 0.15, 0.090, 0.055, and 0.034 of the peak 
pixel. The two lowest contours arc not drawn for the HST, Gemini, VLT (2006), and Keck J- and i'-band images . 



This library served to simulate the unsaturated images 
of LHS 2397aA and LHS 2397aB. To simulate the sat- 
urated images of LHS 2397aA we used other targets 
from the same program that had also been intention- 
ally saturated to search for faint companions. Wc used 
the five targets closest in apparent flux to LHS 2397aA 
(LHS 2924, BRI 0021-0214, TVLM 513-46546, LP 412- 
31, and LHS 2243), which were < 30% different in flux as 
measured by our PSF-fitting routine. We did not scale 
these images to match LHS 2397aA's flux since it is im- 
possible to scale the saturated pixels. Wc added images 
of single stars from our library to these saturated im- 
ages after scaling them down and adding Poisson noise 
to match the S/N of LHS 2397aB. We shifted these im- 
ages by an integer number of pixels that best reproduced 
the observed separation of (+4.1, +4.4) pixels. By apply- 
ing our PSF-fitting routine to these images, we were able 
to assess the position and flux measurements of a faint 
companion made in the presence of a saturated primary. 

From these Monte Carlo simulations, we derived root- 
mean-square (rms) errors in measured {x, y) positions of 
(0.05, 0.06) pix and (0.07, 0.07) pix for the 300 s and 
1000 s images, respectively, with no significant system- 
atic offsets.^ For the fluxes, the ratio of the input to best- 
fit value was 1.05+0.05 for the 300 s image and 1.07+0.11 
for the 1000 s images. In addition, we were able to as- 

Wc found significant statistical correlation between the truth- 
minus-fittcd offsets of the x and y positions {a'^y / a^iyy ~ 0.6). We 
accounted for this covariance when propagating errors to compute 
the uncertainty in the separation and P.A. 



sess the intrinsic uncertainty in our TinyTim PSF-fitting 
routine (i.e., the error when not in the presence of a sat- 
urated star) by comparing each undithered pair of PSF 
library images taken in the same orbit. From this en- 
semble of single-star position and flux measurements we 
found that the scatter in both x and y position measure- 
ments was 0.03 pix, and the scatter in peak flux was a 
factor of 0.06. We adopt these errors for the position 
and flux of LHS 2397aA measured in the 2 s image, and 
we adopt the errors from the Monte Carlo for the posi- 
tions and fluxes of LHS 2397aB measured in the longer, 
saturated-primary exposures. This results in a separa- 
tion uncertainty of 0.08 pix (4 mas), a P.A. uncertainty 
of 0.8°, and a flux ratio uncertainty of 0.08 mag. 

The drift correction we derived is thus quite significant 
compared to our derived positional uncertainties. How- 
ever, if we do not apply this drift correction, the scatter 
in the measured position of LHS 2397aB is twice as large 
as predicted from our Monte Carlo simulations. As dis- 
cussed in Section 3.1, the orbit determination results in 
an almost identical dynamical mass if this single HST 
epoch is excluded. Thus, our derived relative astrometry 
for the HST data is consistent with the orbit as deter- 
mined from all other available astrometry. 

2.2. Keck/NIRC2 LGS 

We have monitored LHS 2397aAB using the laser guide 

star adaptive optics (LGS AO) system at the Keck II 
Telescope on Mauna Kea, Hawaii (Wizinowich et al. 
2006; van Dam et al. 2006). We used the facility near- 
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infrared camera NIRC2 in its narrow field-of-view mode 
because this provides the finest pixel scale. At each 
epoch, we obtained data in one or more filters covering 
the standard atmospheric windows from the Mauna Kea 
Observatories (MKO) filter consortium (Simons & Toku- 
naga 2002; Tokunaga et al. 2002). We also obtained data 
with the CH/^s filter, positioned shortward of a methane 
absorption band (Ac = 1.592 ^m, AA = 0.126 y^m), to en- 
able more robust photometric discrimination in the spec- 
tral type estimate of LHS 2397aB. 

The LGS provided the wavefront reference source for 
AO correction, with the exception of tip-tilt motion. The 
LGS brightness, as measured by the flux incident on the 
AO wavefront sensor, was equivalent to a V ~ 10.4- 
10.7 mag star. The tip-tilt correction and quasi-static 
changes in the image of the LGS as seen by the wave- 
front sensor were measiircd contemporaneously by a 
second, lower bandwidth wavefront sensor monitoring 
LHS 2397a, which saw the equivalent of a i? ~ 15.7- 
16.2 mag star.^ 

On each observing run, we obtained dithered images, 
offsetting the telescope by a few arcscconds between ev- 
ery 2-5 images. The sodium laser guide star was posi- 
tioned at the center of the NIRC2 field-of-view for all 
observations. The images were reduced in a standard 
fashion. We constructed flat fields from the differences 
of images of the telescope dome interior with and with- 
out continuum lamp illumination. Then we created a 
master sky frame and subtracted it from the individ- 
ual images. For L'-band images, the sky subtraction was 
done by pairwise subtraction of consecutive images. Sky- 
subtracted images were registered and stacked to form a 
final mosaic, though all the results described here were 
based on analysis of the individual images. Outlier im- 
ages with much poorer full-width half-maxima (FWHM) 
and/or Strehl ratios were excluded from the analysis. We 
used the publicly available IDL routine NIRC2STREHL'^ to 
estimate the Strehl and FWHM at each epoch, and in 
Table 1 the mean and standard deviation of these values 
over each image set is reported. The variation in Strehl 
ratio and FWHM from between epochs is due to different 
seeing conditions. 

To determine the relative positions and fluxes of 
LHS 2397aA and LHS 2397aB in the imaging data, we 
used a simple analytic representation of the PSF to de- 
blend the two components. The model was the sum of 
three elliptical Gaussians in which each Gaussian compo- 
nent was allowed to have a different FWHM and normal- 
ization, but all components had the same cllipticity and 
semimajor axis P.A. The best-fit parameters were found 
by a Levenberg-Marquardt least-squares minimization in 
which all pixels were weighted equally. This fitting pro- 
cedure yielded a set of measurements of the projected 
separation, P.A., and flux ratio for LHS 2397aAB. We 
used the astrometric calibration from Ghez et al. (2008), 
with a pixel scale of 9.963±0.005 mas/pix and an orien- 
tation for the detector's +y-axis of -|-0.13±0.02 east of 

^ The Keck LGS AO tip-tilt sensor is very red-sensitive, so a red 
source like LHS 2397a appears much brighter than would otherwise 
be expected (R = 17.9 mag, I = 14.6 mag; Monet et al. 2003). For 
late-M dwarfs, we have found that the Keck tip-tilt sensor generally 
sees the equivalent of an (/ -|- 1.0) mag star. 

® http : //www2 . keck . hawaii . edu/optics/lgsao/ software/ 
nirc2strehl . pro 



TABLE 1 
Keck LGS AO Observations 



Date 


Time 


Airmass 


Filter 


FWHM 


Strehl ratio 


(UT) 


(UT) 






(mas) 




2007 Apr 22 


07:46 


1.19 


Ks 


53.8±0.2 


0.277±0.006 




07:56 


1.19 










10:00 


1.47 


Ks' 






2008 Jan 15 


14:48 


1.21 


Ks' 








14:59 


1.23 


Ks 


61.0±0.8 


0.201±0.017 


2009 Jan 23 


12:56 


1.21 


Ks 


71.2±1.1 


0.152±0.015 




13:03 


1.20 


H 


71.9±2.7 


0.059±0.006 




13:14 


1.20 


J 


78.4±4.7 


0.021±0.004 




13:33 


1.19 


L' 


93.4±a.9 


0.434±0.069 



Aperture masking observations. 



north. We applied the distortion correction developed 
by B. Cameron (private communication, 2007) to the as- 
trometry, which changed the results well below the Ict 
level. We also computed the change in the relative as- 
trometry due to differential chromatic refraction (DCR) 
in the same manner as (Dupuy et al. 2009) but did not 
apply this correction to our measurements because we 
found that it induces only a < O.lcr change. 

To assess systematic errors in our PSF-fitting proce- 
dure, we also applied it to simulated Keck images of 
LHS 2397aAB. We scaled down each individual image 
until the primary flux matched the flux from the faint 
companion. In this scaled image, the companion flux is 
at the level of the noise. We then shifted and added this 
image to the original image, matching the binary sepa- 
ration of LHS 2397aAB but avoiding P.A.s within ±60° 
of LHS 2397aB, which were masked out in the original 
image. After running our PSF-fitting routine on these 
images, the resulting scatter in the truth-minus-fitted pa- 
rameters was comparable to or somewhat smaller than 
the rms scatter of the individual measurements (< 0.5c7 
different), thus indicating that the rms scatter is a rea- 
sonable assessment of the errors. The Monte Carlo simu- 
lations also typically indicated significant systematic off- 
sets (> Icr) in the astrometry and also the flux ratio 
(except for the 2009 data). However, we did not ap- 
ply these offsets to any of the parameters derived from 
PSF fitting as these offsets have an insignificant impact 
on the resulting mass estimate (see Section 3.1). They 
also bring the if^-band flux ratios derived from Keck 
imaging out of agreement with the ensemble of flux ratio 
measurements. This is not imexpected since the simu- 
lated images cannot exactly reproduce the science data, 
which has a somewhat asymmetric PSF and speckles that 
cause distinct systematic offsets at different separations 
and P.A.s. 

In 2007 and 2008, data were also obtained using the 
9-hole non-redundant aperture mask installed in the fil- 
ter wheel of NIRC2 (Tuthifl et al. 2006). The data were 
taken in two dither positions separated by 3'.'5, with four 
to six 50 s exposures taken at each dither position. With 
an integrated K-band brightness of 10.7 mag, LHS 2397a 
was fainter than typical targets previously observed with 
adaptive optics and the aperture mask, but it was still 
bright enough not to be limited by readout or background 
noise. Typical interferograms are shown in Figure 2. 
The pipeline used to reduce the aperture masking data 
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was similar to that used in previous papers containing 
NIRC2 masking data (Ireland ct al. 2008; Kraus et al. 
2008; Ireland & Kraus 2008), except that no compara- 
ble single star was observed for calibration of the closure 
phases or squared visibilities. For this reason, we chose 
to only fit to the closure phases, which from previous ex- 
perience with NIRC2 never show an rms scatter of more 
than 3° on bright calibrators in H- or iiT-band. Use of 
squared visibilities would have required a model of the 
fringe decorrelation due to an imperfect AO system, a 
complexity we chose not to tackle. The resulting closure 
phases from our observations arc shown in Figure 2. The 
closure phases arc intrinsically a 4-dimensional function 
(the phase of the bispectrum, e.g., Lohmann et al. 1983), 
so they are difficult to present on a 2- dimensional plot. 
Note that the closure phases are a linear function of con- 
trast for contrast ratios larger than ^5:1, so for example 
if the companion were twice as bright, then all the mea- 
sured and model closure phases would be doubled. 

The closure phase uncertainties were initially approxi- 
mated by the standard error of the mean calculated from 
the scatter among individual exposures. The uncertain- 
ties were subsequently increased by adding a calibration 
error in quadrature so that the resulting reduced of 
the fit was 1.0. Although we fit all 84 closure phases 
from the 9-hole mask, only 28 of these are formally in- 
dependent. To correctly account for this non-diagonal 
covariance matrix in our binary fitting, we scaled the er- 
rors in the least-squares fit to the data by ^84/28. This 
process has been validated both by a comparison to fits 
using full covariance matrices (Kraus et al. 2008) and by 
orbit fits using mixed data that resulted in a reduced 
consistent with unity, where the orbit fit had many 
degrees of freedom (e.g., Martinache et al. 2007). 

In Figure 2, it can be seen that binary is clearly de- 
tected at > 3a in many individual closure triangles for 
the two Ks-hand detections. The C-ff4.s-band detection 
has only one closure triangle above 3a but nonetheless 
has a very clear and unique fit to all triangles taken to- 
gether that is consistent with the other data from the 
same epoch. 

The astromctry and flux ratios derived from the Keck 
data are given in Table 2. The different data sets at 
each epoch, both aperture masking and direct imaging, 
give consistent measurements of the binary parameters. 
When fitting the orbit, we adopt the measurements from 
the single data set at each epoch with the smallest astro- 
metric errors. For the 2008 data, the error in separation 
is smaller for the masking data while the P.A. error is 
smaller in the imaging data. Because the binary separa- 
tion is very close to the outer limit at which masking ob- 
servations are unambiguous in the 2008 data (wl50 mas; 
imposed by the smallest non-redundant baselines in the 
mask) , we conservatively chose to adopt imaging astrom- 
etry for the 2008 epoch when fitting the orbit. Regard- 
less of which Ks-ha.nd data set we use at each epoch, the 
resulting mass estimate is not significantly changed, as 
discussed in Section 3.1. 

2.3. Gemini/Hokupa'a 

LHS 2397aAB was imaged on UT 2002 February 7 

by the Hokupa'a curvature AO system at the Gemini- 
North Telescope on Mauna Kea, Hawai'i. We retrieved 
these raw data from the Gemini science archive and reg- 
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Fig. 2.— Keck interferograms of LHS 2397aAB ob- 
tained using NIRC2's 9-hole aperture mask (square-root 
stretch). The left panels show the modeled versus mea- 
sured closure phases (C.P.) at each epoch. The bi- 
nary parameters were derived from the modeled closure 
phases, and their errors were assessed in a Monte Carlo 
fashion that accounted for the measured closure phase 
errors. 

istered, sky-subtracted, and performed cosmic ray rejec- 
tion on the images. Figure 1 shows a typical image from 
of one of the 15 JiT'-band 3 s exposures which were used 
to derive the astrometry for LHS 2397aAB. Analysis of 
these data has previously been presented by Freed et al. 
(2003); however, astrometric errors were not derived in 
that work, so we conducted our own analysis of these 
data. 

We used the same analytic PSF-fitting routine as 

for the Keck data to fit the Gemini images of 
LHS 2397aAB. Adopting the nominal instrument pix- 
elscale of 19.98±0.08 mas/pix, we found a separation of 
205±8 mas, where the uncertainty is the standard de- 
viation of measurements from individual dithers. This 
is in good agreement with the 207 mas separation re- 
ported by Freed et al. (2003). However, the P.A. we 
find (206!8±2.4°) is significantly different from the 152.0° 
P.A. reported by Freed et al. (2003). This must be due to 
an inaccurate reporting of the orientation of the camera 
in the header of the archival data, which indicates that 
the P.A. of the +y-axis is 0° . Such a discrepancy was pre- 
viously seen in our work with archival Gemini/Hokupa'a 
data of HD 130948BC and Gl 569Bab (Dupuy et al. 
2009). Thus, we adopt the published value of the P.A., 
but not its quoted error, for our analysis. 

To assess systematic errors in the Gemini astrometry, 
we simulated many Gemini images of LHS 2397aAB in 
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TABLE 2 

Best-Fit Binary Parameters for LHS 2397aAB 



1--, 1 / Y Tm \ 

Epoch (UT) 


Instrument 


Filter 


p (mas) 


P.A. (°) 


Am (mag) 


1997 Apr 12 


JTST/WFPC2-PC1* 


F814W 


274±4 


87.3±0.8 


4.18±0.08 




Gemini /Hokupa' a 


i\ 


ZUOzto 


1 CO 0-1-0 A 


O.UDztU.40 


2003 May 31 


VLi /N AGO 
VLi /IN AUU 


Ks 


168±8 


188.D±1.2 


2.75±0.16 


2005 Jan 12 


Ks 


129±5 


27d.2±1.4 


2.84±U.lD 


2007 Apr 22 


Keck/NIRC2 


Ks 


119.7±2.9 


348.4±0.8 


2.94±0.10 




Keck/NIRC2 (masking) 


CHas 


112.5±2.8 


347.0±1.3 


3.37±0.24 




Keck/NIRC2 (masking)* 


Ks 


115.3±1.8 


349.3±1.1 


2.78±0.11 


2008 Jan 15 


Keck/NIRC2 (masking) 


Ks 


144.6±2.3 


24.1±0.7 


2.97±0.11 




Keck/NIRC2* 


Ks 


146±4 


25.1±0.5 


2.70±0.12 


2009 Jan 23 


Keck/NIRC2* 


Ks 


204.0±1.5 


53.8±0.3 


2.80±0.03 




Keck/NIRC2 


H 


202.0±2.4 


53.9±0.5 


2.96±0.05 




Keck/NIRC2 


J 


196±5 


53.9±0.8 


3.12±0.08 




Keck/NIRC2 


L' 


200±3 


53.3±0.4 


1.92±0.06 



*Used in the orbit fit. 



the same fashion as the Keck data and fit each one with 
our analytic PSF model. The resulting scatter in the 
truth-minus-fitted parameters was 14 mas in separation 
and 4° in P.A., with no significant systematic offsets. 
These errors are much larger than the rms scatter of the 
individual measurements, which is likely because the PSF 
in the science data is somewhat asymmetric so that the 
offset between the best-fit and input parameters varies 
widely, depending on the binary P.A., and by averag- 
ing over many binary P.A.s we overestimate the scatter. 
However, as shown in Section 3.1, whether we adopt the 
rms scatter or the Monte Carlo errors only affects the 
of the orbit and not the resulting mass estimate. In fact, 
even if the Gemini epoch is excluded from the orbit fit 
the resulting mass estimate changes insignificantly. 

2.4. VLT/NACO 

We retrieved archival images of LHS 2397aAB ob- 
tained with the Very Large Telescope (VLT) at Paranal 
Observatory on UT 2003 May 31 and 2006 Jan 15. These 
data were taken with the NACO adaptive optics sys- 
tem (Lenzen et al. 2003; Rousset et al. 2003) using the 
N90C10 dichroic at both epochs and the S27 and S13 
cameras in 2003 and 2006, respectively. The nominal 
pixel scales of these cameras are 27.053±0.019 mas/pix 
and 13.221±0.017 mas/pix.'^ We registered, sky- 
subtracted, and performed cosmic ray rejection on the 
raw archival images. The 2003 data comprise five Ks- 
band 5 s exposures, and the 2006 data comprise 32 Ks- 
band 10 s exposures. Typical images from each data set 
are shown in Figure 1. 

We used the same analytic PSF-fitting routine as for 
the Keck data to fit the VLT images of LHS 2397aAB. 
From the 2003 VLT data we derived a separation of 
168±8 mas, a P.A. of 188.6±1.2°, and a ifg-band flux 
ratio of 2.75±0.16 mag. From the 2006 VLT data we de- 
rived a separation of 129±5 mas, a P.A. of 276.2±1.4°, 
and a Ks-hand flux ratio of 2.84±0.16 mag. The quoted 
uncertainty is the standard deviation of measurements 

http : //www . eso . org/sci/f acilities/p5iranal/instrumeiit s/ 
naco/doc/VLT-MAN-ESD-14200-2761_v83 . 3 . pdf 



from individual dithers. To assess any additional system- 
atic errors in the VLT astrometry, we simulated many 
VLT images of LHS 2397aAB in the same fashion as the 
Gemini and Keck data. These simulations yielded equiv- 
alent errors to the rms scatter quoted above (< 0.2a 
different), and thus the quoted errors are a good repre- 
sentation of the uncertainties. 

3. RESULTS 

3.1. Orbit Determination & Dynamical Mass 

The orbit of LHS 2397aAB is very well constrained as 
our observations cover most of the orbital period. In 
order to search for the influence of parameter degenera- 
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Fig. 3. — Relative astrometry for LHS 2397aAB along 
with the best fitting orbit (reduced of 1.04 for 7 de- 
grees of freedom). The empty circles show the predicted 
location of LHS 2397aB in the future. Error bars are 
comparable to or smaller than the plotting symbols. The 
orbit is very well constrained, so the uncertainty in the 
total mass is dominated by the 3.0% parallax error. 
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Fig. 4. — Probability distributions of all orbital parameters derived from the MCMC analysis: semimajor axis (a), 
orbital period (P), eccentricity (e), inclination (?'), epoch of periastron (To), P. A. of the ascending node (fi), and 
argument of periastron (lo). Each histogram is shaded to indicate the 68.3% and 95.4% confidence regions, which 
correspond to la and 2a for a normal distribution, and the solid vertical lines represent the median values. Note that 
To is shown in days since UT 2006 Apr 19 12:00 for clarity. 



cies in our orbit fit and determine robust confidence lim- 
its on the orbital parameters, we used a Markov Chain 
Monte Carlo (MCMC) technique (e.g., Bremaud 1999) 
for orbit fitting, in addition to a gradient descent tech- 
nique. In short, the MCMC method constructs a se- 
ries of steps through the model parameter space such 
that the resulting set of values (the "chain") is asymp- 
totically equivalent to the posterior probability distribu- 
tion of the parameters being sought. The code that per- 
formed the MCMC fit is described in detail in the study 
of 2MASS J1534-2952AB by Liu et al. (2008). Chains 
all had lengths of 2x 10* steps, and the correlation length 
of our most correlated chain, as defined by Tegmark et al. 
(2004), was 230 for the argmnent of periastron. This 
gives an effective length of the chain of 8.7x10'"', which 
in turn gives statistical uncertainties in the parameter 
errors of about l/\/8.7 x 10^ = 0.11%, i.e., negligible. 

We used uniform priors in period (P), semimajor axis 
(a), P. A. of the ascending node (il), argument of perias- 
tron (cij), and time of periastron passage (Tg). We used 
a prior in inclination proportional to sin(z) (i.e., random 
orbital orientation) and an eccentricity prior of /(e) = 2e 
(e.g., see Duquennoy & Mayor 1991). Figure 4 shows the 
resulting MCMC probability distributions for the seven 
orbital parameters of LHS 2397aAB. The best-fit param- 
eters and their confidence limits are given in Table 3, and 
the best-fit orbit is shown in Figure 3. The reduced 
of the orbital solution is 1.04, with 7 degrees of freedom. 

Applying Kepler's Third Law to the period and semi- 
major axis distributions gives the posterior probability 



distribution for the total mass of LHS 2397aAB, with a 
median of 0.146 Mq, a standard deviation of 0.004 Mq, 
and 68.3(95.4)% confidence hmits of to.oMl^o.oos) M© 
(Figure 5). The MCMC probability distribution of the 
total mass does not include the uncertainty in the par- 
allax (3.0%), which by simple propagation of errors 
would contribute an additional 9.0% uncertainty in mass. 
In fact, since the MCMC-derived mass distribution is 
slightly asymmetric, we account for this additional er- 
ror by randomly drawing a normally distributed paral- 
lax value for each step in the chain, which we then used 
to compute the total mass. The resulting mass distri- 
bution is nearly indistinguishable from Gaussian (Fig- 
ure 5). Our final determination of the total mass is 
146+°°i''' 

U.i'iD_o.oi3 



(-0:024) M© at 68.3(95.4)% confidence. 



As an independent verification of our MCMC results, 
we also fit the orbit of LHS 2397aAB using the linearized 
least-squares routine ORBIT (described in Forveille et al. 
1999). All of the orbital parameters and their errors are 
consistent between the ORBIT and MCMC results, as ex- 
pected for well constrained parameters with nearly Gaus- 
sian probability distributions. The and total mass of 
the ORBIT solution and the MCMC solution were identi- 
cal. 

From the visual orbit alone there is a 180° ambigu- 
ity in the P. A. of the ascending node (O), and corre- 
spondingly Lu, that can only be resolved by radial velocity 
measurements. Basri & Reiners (2006) measured the ra- 
dial velocity of LHS 2397a at two epochs (UT 1995 Mar 
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TABLE 3 

Derived Orbital Parameters for LHS 2397aAB 



MCMC ORBITt 



Parameter 


Median 


68.3% c.l. 


95.4% c.l. 




Semimajor axis a (mas) 


216.4 


-1.9, 2.0 


-3.8, 3.9 


216.5±2.0 


Orbital period P (yr) 


14.22 


-0.10, 0.10 


-0.19, 0.20 


14.22±0.10 


Eccentricity e 


0.350 


-0.005, 0.005 


-0.011, 0.011 


0.345±0.005 


Inclination i (°) 


41.9 


-1.1, 1.1 


-2.3, 2.2 


42. Oil. 1 


Time of poriastron passage To - 2453845.0^ (JD) 





-15, 14 


-30, 28 


0±15 


P.A. of the ascending node Q. (°) 


257.8 


-1.5, 1.5 


-3.1, 2.9 


257.8±1.5 


Argument of periastron w (°) 


38.5 


-2.3, 2.4 


-4.6, 4.8 


38.5±2.4 


Total mass (Mq): fitted'" 


0.146 


-0.004, 0.004 


-0.008, 0.008 


0.146±0.004 


Total mass (M©): finaF 

Reduced (7 degrees of freedom) 


0.146 


-0.013, 0.015 


-0.024, 0.031 


0.146±0.14 


1.04 






1.04 



tThe orbital parameters determined by ORBITwith their linearized la errors. 
'^UT 2006 Apr 19 12:00:00.0 

''The "fitted" total mass represents the direct results from fitting the observed orbital motion of the two 
components without accounting for the parallax error. For the linearized ORBIT error, the covariance between 
P and a is talcen into account. 

''The "final" total mass includes the additional 9.0% error in the mass due to the error in the parallax. 
This final mass distribution is essentially Gaussian. 



12 and 2002 May 19) in the optical where LHS 2397aB 
is essentially invisible. For a mass ratio of 0.7 (Sec- 
tion 4.2), the expected difference in the radial velocity of 
the primary between these two epochs is ±1.8 km s~^, 
where the unknown sign reflects the 180° ambiguity in CI. 
Basri & Reiners (2006) measured a velocity difference of 



T — ' — ' — ' — r 




0.12 0.14 0.16 0.18 
Totol moss (Mq) 

Fig. 5. — Probability distribution of the total mass of 
LHS 2397aAB resulting from our MCMC analysis. The 
histogram is shaded to indicate the 68.3%, 95.4%, and 
99.7% confidence regions, which correspond to la, 2a, 
and 3a for a normal distribution. The dashed line rep- 
resents the median value of 0.146 Mq. The standard 
deviation of the distribution is 0.004 M0. The dotted 
unshaded curve shows the final mass distribution after 
accounting for the additional 9.0% error due to the uncer- 
tainty in the parallax; the result is essentially Gaussian. 
The confidence limits for both distributions are given in 
Table 3. 



+ 1.1±1.6 km s^^, which is consistent with +1.8 km s"-*^ 
but discrepant with —1.8 km at 1.8<t. 

We have tested whether varying the input astrom- 
etry and corresponding uncertainties affects the orbit- 
fitting results using ORBIT. Regardless of what astrom- 
etry was used, the mass estimate and its uncertainty 
did not change significantly. We tested many permu- 
tations using ORBIT, such as using Keck imaging instead 
of masking astrometry or Monte Carlo instead of rms 
errors, and we found that the resulting mass estimates 
ranged from 0.142-0.150 Mq and uncertainties ranged 
from 0.003-0.006 Mq. The of these different permu- 
tations ranged from 5-17 (7 degrees of freedom), so many 
of these scenarios are clearly not as favorable as our de- 
fault solution that has a x^ of 7.27. The permutations 
with lower x^ simply correspond to the cases with some- 
what larger assumed astrometric uncertainties, and these 
had essentially identical masses and mass errors as our 
default solution. 

In fact, any single epoch can be excluded and the re- 
sulting dynamical mass changes by < 0.006 Mq, with 
the uncertainty from the orbit alone increasing to only 
±0.006 Mq, which is well below the error due to the par- 
allax uncertainty. Furthermore, pairs of epochs can even 
be excluded without significantly changing the dynam- 
ical mass. For example, both of the HST and Gemini 
epochs or both of the VLT epochs can be excluded and 
the resulting dynamical mass changes by < 0.008 Mq 
and the error from the orbit alone increases to 0.013 Mq 
(i.e., comparable to the error due to the parallax; uncer- 
tainty). Thus, our mass estimate for LHS 2397aAB and 
its uncertainty is insensitive to the input astrometry and 
the method of determining the measurement errors. 

3.2. Spectral Types 

Spectral type determination for the components of 
LHS 2397aAB is impeded by the lack of resolved spec- 
troscopy. However, the optical classification of M8 (Kirk- 
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Fig. 6. — Measurements of the projected separation (lejf) and P.A. (right) of LHS 2397aAB. The best-fit orbit is 
shown as a sohd Hne. The bottom panels show the observed minus predicted measurements with observational error 
bars. The fact that the 2007 Keck astrometry shows the smallest errors even though the data were obtained near the 
binary's minimum separation highlights the exquisite precision that can be obtained for such a tight, high-contrast 
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Patrick et al. 1995) can be readily applied to the pri- 
mary since the /-band flux ratio of « 4 mag indicates 
that LHS 2397aB contributes negligibly to the flux in 
the unresolved optical spectrum. The spectral type of 
LHS 2397aB must be estimated indirectly, and we have 
done so in a variety of ways. 

The available IJHKL' colors of LHS 2397aB do not 
provide a strong constraint on the spectral type; however, 
they are consistent with late-L dwarfs (Dahn et al. 2002; 
Golimowski et al. 2004; Liebert & Gizis 2006). Unfortu- 
nately, even if the colors of LHS 2397aB were measured 
perfectly, the dispersion among known L dwarfs limits a 
spectral type estimate from colors alone to at least ±2 
subclasses. 

We also used the empirical absolute magnitude re- 
lations of Liu et al. (2006)* with known and sus- 
pected binaries removed to estimate the spectral type 
of LHS 2397aB. We converted these relations from the 
MKO to the 2MASS photometric system using the re- 
lations of Stephens & Leggett (2004). In addition, we 
derived the L'-band absolute magnitude relation for the 
same sample, though 7 of the 29 objects from Liu et al. 
(2006) could not be used as they did not have published 
L'-band photometry (Leggett et al. 2002; Golimowski 
et al. 2004). The resulting polynomial fit was 

Mi' =9.759 + 0.1507 x SpT + 2.868 x 10"^ x SpT^ 

- 3.674 X 10-^ X SpT^ + 1.338 x IQ-'^ x SpT'' 

with an rms scatter about this relation of 0.18 mag. 
By simply interpolating these absolute magnitude re- 
lations, the if-band photometry of LHS 2397aB gave a 
spectral type of L7.5±0.3, the 77-band photometry gave 
L6.6±0.3, and the L'-band photometry gave L7.4±0.5. 
These quoted uncertainties only account for the photo- 
metric errors (i.e., the la spectral type range corresponds 
to the 1(7 range in photometry) and not scatter in the 
empirical relations. Because of the flattening of the ab- 
solute magnitude relation in J-band, we were unable to 

* These empirical relations are based on near-infrared spectral 
types. 



estimate a robust spectral type from the J-band pho- 
tometry, although it is consistent with the spectral type 
ranges listed above. 




L5 TO T5 
Spectral Type 

Fig. 7. — The final probability distribution of the spec- 
tral type of LHS 2397aB (solid black hne), derived using 
the absolute magnitude relations of Liu et al. (2006), with 
suspected binaries removed, and our own L'-band rela- 
tion as described in the text. The joint probability dis- 
tribution of the spectral type was derived from J-band, 
iJ-band, X-band, and L'-band photometry (thin colored 
lines). The filled circle with error bars denotes the me- 
dian and 1(7 confidence limits of this final spectral type 
distribution. 

We also employed a more sophisticated method for 
spectral type estimation that accounts for both the un- 
certainty in the absolute magnitude of LHS 2397aB and 
the scatter in the empirical relations. At each half-integer 
spectral type, we considered two Gaussian probability 
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TABLE 4 

Measured Properties of LHS 2397aAB 



Jr^rop6rty 


LHS 2397aA 


LHS 2397aB 


Ref 


Mtot (Mq) 


0.146±0.014 


1 


Semimajor axis (AU) 


3.09±0.10 


1,2 


d (pc) 


14.3± 0.4 


2 


Spectral Type 


M8±0.5 


L7±l 


1,3 


J (mag) 


12.11±0.02 


15.23±0.08 


1,4 


H (mag) 


11.26±0.02 


14.22±0.05 


1,4 


K (mag) 
L' (mag) 


10.83±0.02 


13.64±0.04 


1,4 


10.20±0.02 


12.12±0.05 


1,5 


J — H (mag) 


0.85±0.03 


l.OliO.lO 


1,4 


H-K (mag) 


0.43±0.03 


0.58±0.07 


1,4 


J — K (mag) 
K-L' (mag) 


1.28±0.03 


1.59±0.09 


1,4 


0.63±0.03 


1.51±0.07 


1,4,5 


Mv (mag) 


18.80±0.07 




2 


Mj (mag) 


11.34±0.07 


14.46±0.10 


1,2,4 


Mh (mag) 


10.49±0.07 


13.45±0.08 


1,2,4 


Mk (mag) 


10.06±0.07 


12.87±0.08 


1,2,4 


Ml, (mag) 
log(Lboi/L0) 


9.43±0.07 


11.35±0.09 


1,2,5 


-3.42±0.03 


-4.52±0.03 


1 


Alog(Lboi) 


1.102±0.011 


1 



References. — (1) This work; (2) Monet et al. (1992); 
(3) Kirkpatrick et al. (1995); (4) Cutri et al. (2003); 
(5) Leggett et al. (2002). 

^AU near-infrared photometry on the MKO system. 



distributions: (1) the distribution of the observed abso- 
lute magnitude of LHS 2397aB and (2) the distribution 
of absolute magnitudes predicted by the empirical rela- 
tion from the overall scatter in the relation as given by 
Liu et al. (2006). We multiplied these two distributions 
and integrated the result in order to obtain a measure 
of the relative probability of LHS 2397aB having that 
spectral type. If the two distributions overlapped sub- 
stantially, then the integral was large, but if the two dis- 
tributions were very discrepant the integral went to zero. 
Thus, we built an ensemble of relative probabilities which 
we adopted as the probability distribution of the spectral 
type of LHS 2397aB (Figure 7). We did this separately 
for each bandpass, and the derived spectral type and un 
certainty was T0.9t4;g for -'^-band, L7.2 



+i;3 for 7J-band, 



L7.6li;| for if-band, and L7.4+};i for L'-band (68.3% 
confidence limits) . Since each bandpass represents an in- 
dependent constraint on the spectral type, we took the 
product of these different spectral type distributions to 
determine the final spectral type probability distribution 
(Figure 7). This is conceptually equivalent to taking the 
weighted mean of the spectral types listed above, and 
the result is L7.2±0.7.^ Because of the good agreement 
between the spectral types inferred from different band- 
passes, the final derived spectral type of LHS 2397aB is 
essentially the same with smaller errors. 

We adopt the best spectral type estimates available: 
for LHS 2397aA this is M8±0.5 from the optical spec- 
trum, and for LHS 2397aB we assign a spectral type of 
L7±l from the absolute magnitude relations, where the 

^ If we use the Liu et al. (2006) empirical relations which 
do not exclude suspected binaries, the resulting spectral type of 
LHS 2397aB is L7.6lJ;g. 
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Fig. 8. — The integrated light spectrum of LHS 2397a, 
with optical data (< 0.9 ^m) taken from Kirkpatrick 
et al. (1995) and our near-infrared data obtained with 
IRTF/SpeX. Data for M dwarf spectral standards are 
shown for comparison (Henry & Kirkpatrick 1990; Bur- 
gasser et al. 2004; Burgasser & McElwain 2006; Bur- 
gasser et al. 2008). LHS 2397a is optically typed as M8, 
and its near-infrared {H and K band) spectrum is dis- 
crepant with the standard VB 10 due to additional flux 
from the very red L7 dwarf companion LHS 2397aB. 

uncertainty accounts for both the observational errors in 
the photometry and the scatter in the empirical relations. 

3.3. Bolometric Luminosities 

We computed bolometric luminosities for 
LHS 2397aAB in two different ways: (1) directly 
using integrated-light measurements of LHS 2397a and 
(2) using our resolved photometry for the individual 
binary components along with spectral templates. 
Kirkpatrick et al. (1995) measured the optical integrated 
-light spectrum of LHS 2397a («0. 64-0.92 nm)}° On 
2008 June 27 UT we obtained a near-infrared spec- 
trum of LHS 2397a («0.82-2.4 ^m) using IRTF/SpeX 
(Rayner et al. 1998) in SXD mode (i? = 1200), which 
we reduced using the SpeXtool software package (Vacca 
et al. 2003; Gushing et al. 2004). We combined these 
spectra, using optical data below 0.9 fim and our 
SpeX data above 0.9 fim. Figure 8 shows the resulting 
integrated-light spectrum of LHS 2397a, along with 
M dwarf spectral standards from Kirkpatrick et al. 
(1991)." 

We derived a total bolometric luminosity of 
log(Lboi/L0) = -3.33±0.03 for LHS 2397a using 
the integrated-light spectrum (0.6-2.4 fim), L'-band 
(3.8 /im) photometry from Leggett et al. (2002), and 

Publicly available at http://www.dwarfarchives.org. 

We combined optical and near-infrared data for the M dwarf 
spectral standards in the same way as LHS 2397a. The optical data 
are from Henry & Kirkpatrick (1990), and the near-infrared data 
are SpeX prism spectra from Burgasser et al. (2004), Burgasser & 
McElwain (2006), and Burgasser et al. (2008). 
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Fig. 9. — The heliocentric space velocity of LHS 2397a (star) is shown alongside other ultracool dwarfs: > M7 dwarfs 
(squares), L dwarfs (circles), and T dwarfs (triangles). The la and 2(7 ellipsoids of the thin disk (dashed) and thick 
disk (dotted) as predicted by the Besangon galaxy model (Robin et al. 2003) are also shown for comparison. The 
space velocity of LHS 2397a is consistent with other ultracool dwarfs, and we derive a 99.2% thin disk membership 
probability. 



24 /im Spitzer/MIPS photometry from Gautier et al. 
(2007). We neglected any flux at shorter wavelengths, 
interpolated between gaps in the data, and extrapolated 
the flux beyond 24 /im assuming a blackbody. We 
determined the luminosity error in a Monte Carlo 
fashion by adding randomly drawn noise to our data 
over many trials and computing the rms of the resulting 
luminosities. 

For LHS 2397a, the SpeX spectrum and L'-band pho- 
tometry alone account for 98% of its luminosity. Thus, 
it is possible to accurately estimate the luminosities of 
the individal components using only our resolved JHKL' 
photometry and spectral templates to interpolate prop- 
erly between bandpasses. We selected objects from 
the SpeX Prism Libraries^^ that had similar colors to 
LHS 2397aA and LHS 2397aB and calibrated them to 
match our resolved JHK photometry. In fact, we found 
that the exact choice of the templates had an insigniflcant 
impact on the resulting luminosity measurements. We 
integrated over the template spectra, neglecting flux at 
shorter wavelengths, interpolating the flux between the 
spectra and L' band, and extrapolating to longer wave- 
lengths assuming a blackbody. Thus, we derived lumi- 
nosities of log(Lboi/L0) = -3.42±0.03 for LHS 2397aA 
and -4.52±0.03 for LHS 2397aB. The sum of these lumi- 
nosities is — 3.39±0.03, which is consistent with the Lboi 
derived from integrated-light data. The dominant source 
of error in these luminosities is the parallax, which is 
common to both components, and in the following anal- 
ysis we correctly account for this co variance. 

3.4. Age Constraints from Kinematics and Activity 

The space motion of LHS 2397a could provide a con- 
straint on its age, especially if it were found to be kine- 
matically associated with a population of known age 
(e.g., the thick disk). Using the proper motion and 
parallax from Monet et al. (1992) and a radial veloc- 
ity of 34±2 km s'^ (Mohanty & Basri 2003; Tinney & 
Reid 1998), we computed the heliocentric velocity of 

http: / /www . br owndwar f s . org/spexprism/ 



LHS 2397a to be iU,V,W) = (-26.5±0.8, -36.1±1.5, 
+ 12.5±1.4) km s"^."'^^ We adopt the sign convention for 
U that is positive toward the Galactic center and account 
for the errors in the parallax, proper motion, and radial 
velocity in a Monte Carlo fashion. 

For comparison, we compiled all the radial velocities 
for objects with spectral types of M7 or later from Mo- 
hanty & Basri (2003), Bailer- Jones (2004), and Zapa- 
tero Osorio et al. (2007) and computed space velocities 
for those objects with parallaxes and proper motions 
from Dahn et al. (2002), Monet et al. (1992), Tinney 
et al. (2003), Vrba et al. (2004), Tinney (1996), Tinney 
et al. (1995), and Henry et al. (2006). The resuhing he- 
liocentric velocities are shown in Figure 9. The mean 
and rms scatter of the space velocity of this population 
is {U,V,W) = (-18±37, -24±26, -5±15) km s'^, and 
LHS 2397a is only 1.3cr away from the mean of this ellip- 
soid. Thus, the space motion of LHS 2397a is not signifi- 
cantly different from other ultracool dwarfs, implying an 
age consistent with the population of ultracool dwarfs as 
a whole. Several authors have attempted to estimate the 
age of this population, typically comparing the distribu- 
tion of tangential velocities (Vtan, which requires only a 
proper motion and distance determination) to the well 
studied nearby populations of FGKM stars. The result- 
ing age for the population of ultracool dwarfs estimated 
in this way has been found to be 2-4 Gyr (Dahn et al. 
2002; Faherty et al. 2009)." 

We have also assessed LHS 2397a's membership in the 

This is consistent with the heliocentric velocity derived by 
Freed et al. (2003), but inconsistent with the velocity computed 
by Tinney & Reid (1998). It appears that their space velocity 
for LHS 2397a is discrepant by ^2x{U, V, W)q as a result of sub- 
tracting rather than adding the solar motion to the heliocentric 
velocities. 

^* Zapatero Osorio et al. (2007) determined a somewhat younger 
age (~1 Gyr) for the population of ultracool dwarfs, based on the 
small sample of L and T dwarfs will full space velocities (21 ob- 
jects). However, since L and T dwarfs span a wider range of masses 
than earlier type objects, a typical IMF that rises at lower masses 
will naturally increase the number of young objects in this sample, 
biasing a kinematically derived age (e.g., see Section 4.5 of Dahn 
et al. 2002). 
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Fig. 10. — Total mass (Mtot) predicted by evolutionary models as a function of age, given the observational constraint 
of the luminosities of the two components of LHS 2397aAB. The curved shaded regions show the la and 2(T ranges in 
this model-derived mass. By applying the additional constraint of the measured total mass (Mtot), we used the models 
to determine the age of LHS 2397aAB (see Section 4.1). The horizontal gray bars show our Icr and 2(T constraints 
on the total mass, and the resulting median, Icr, and 2a model-inferred ages are shown by solid, dashed, and dotted 
lines, respectively. Model-inferred ages are truncated at 10 Gyr (the oldest age included in both sets of models) - this 
happens at Icr for the Lyon models and 2cr for the Tucson models. This is because the high mass tail of the Mtot 
distribution corresponds to both LHS 2397aA and LHS 2397aB being stars at the bottom of the main sequence, and 
evolutionary models do not reach faint enough luminosities for such high mass objects by 10 Gyr. 



Galactic populations of the thin disk (1-10 Gyr; e.g., 
Wood & Oswalt 1998) and thick disk (~10 Gyr; e.g., 
Ibukiyama & Arimoto 2002) using the Besangon model of 
the Galaxy (Robin et al. 2003). Because the scale height 
of these structures are much larger than the measured 
distances to ultracool dwarfs, we found that a special- 
ized model population for the volume within 1 square- 
degree on the sky and 5 pc in the radial direction of 
LHS 2397a had essentially the same space motions as 
the model population of all objects within 50 pc over the 
whole sky. Robin et al. (2003) describe the velocity ellip- 
soids of the Besangon model thin and thick disks (shown 
in Figure 9) and how these space motions are allowed 
to evolve over time due to dynamical heating. To deter- 
mine membership probabilities of a given object belong- 
ing to the thin and thick disks, we simply counted the 
nearest 10^ model objects in ([/, V, Vl^)-space and com- 
puted the fraction that belonged to each population. For 
LHS 2397a, we found a membership probability of 99.2% 
for the thin disk and 0.8% for the thick disk.^^ 

Finally, the fact that LHS 2397a is an active Hct flare 
star {\og{Liia I Lho\) ~ —4.22; Gizis ct al. 2000) could also 
potentially provide an age constraint, as the activity of 
M dwarfs changes with age. West et al. (2008) showed 
that the fraction of active M dwarfs as a function the 
vertical distance above the Galactic plane {z) provides a 
constraint on the activity lifetime of M dwarfs, given a 
model of how thick disk heating pumps up z over time. 
West et al. (2008) found that the activity lifetime in- 

We also computed membership probabilities for all the ul- 
tracool dwarfs shown in Figure 9 and found only three systems 
with > 10% probability of belonging to the thick disk: ESO 207-61 
(69%), GJ 1001 (40%), and Gl 802 (36%). 



creases monotonically with M dwarf spectral type, and 
the latest type for which they were able to determine a 
robust hfetime was M7 (8.0t'j';o Gyr). This provides a 
weak constraint on the age of LHS 2397a, as its activity 
is therefore expected to last for at least >8 Gyr. 

In summary, the space motion of LHS 2397a indicates 
that it is (1) consistent with the field ultracool dwarf pop- 
ulation and (2) most likely (99.2%) a member of the thin 
disk. These properties do not strongly constrain its age: 
the age of the ultracool field dwarf population is poorly 
constrained by current observations, and the thin disk is 
believed to have been forming stars continuously over the 
past ~10 Gyr. The activity of LHS 2397a is expected to 
be long-lived, and so it also does not strongly constrain 
the age. Thus, the Icr age range of 1.0-10 Gyr derived 
from evolutionary models for LHS 2397a (Section 4.1) is 
consistent with the broad range of ages allowed by its 
space motion and activity. 

4. TESTS OF MODELS 

Mass is the primary input parameter for evolutionary 
models of ultracool dwarfs. Thus, the strongest tests 
of these models are made possible by objects with a di- 
rect observational constraint on their mass. In the fol- 
lowing analysis, we have used our derived total mass of 
LHS 2397aAB along with its other measured properties 
to test theoretical models. We drew the measured prop- 
erties of LHS 2397aAB from appropriate random dis- 
tributions, carefully accounting for the covariance be- 
tween different quantities (e.g., Mtot and Lboi are corre- 
lated through the distance). We have chosen to use Lboi 
rather than T^s as the basis of our model comparisons 
because values of Toff in the literature are invariably tied 
to either evolutionary or atmospheric theoretical mod- 
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els in some way due to the lack of direct radius mea- 
surements for brown dwarfs. By using iboi, which only 
depends on direct measurements of SEDs and distances, 
we have avoided circular comparisons. This approach fol- 
lows our previous work for 2MASS J1534-2952AB (Liu 
et al. 2008) and HD 130948BC (Dupuy et al. 2009). 

In the following, we consider two independent sets of 
evolutionary models; the Tucson models (Burrows et al. 
1997) and the Lyon DUSTY models (Chabrier et al. 
2000), which are appropriate for both components of 
LHS 2397aAB, as their photospheres are expected to be 
in the effective temperature range over which appreciable 
amounts of dust exist. 

4.1. Model- Inferred Age 

As described in detail by Liu et al. (2008) and Dupuy 
et al. (2009), the total mass of a binary along with its in- 
dividual component luminosities can be used to estimate 
the age of the binary from evolutionary models. This 
age estimate can be surprisingly precise when both com- 
ponents are likely to be substellar since their luminosi- 
ties depend very sensitively on age. For LHS 2397aAB, 
which only has one substellar component, we derive an 
age of l.Slg'e from the Tucson models and an age 
of 1.8^0^ Gyr from the Lyon models (Figure 10). 

Note that there is an upper limit of 10 Gyr on the ages 
estimated from evolutionary models as this is the oldest 
age included in both sets of models. Our age estimates 
for LHS 2397aAB hit this limit at la for the Lyon mod- 
els because the la upper limit on total mass (0.160 M0) 
corresponds to both LHS 2397aA and LHS 2397aB be- 
ing stars at the bottom of the main sequence, and the 
Lyon models do not reach faint enough luminosities for 
such high mass objects by 10 Gyr. This is also the case 
for the Tucson models, though the limit is reached be- 
yond the la upper limit in the total mass. Using models 
that extend to older ages (e.g., the age of the universe) 
would not remedy this problem because LHS 2397aB's 
luminosity is much too low for a star even though the 
total mass formally allows it to be one (see Section 4.2). 

4.2. Individual Masses 

By measuring the relative orbit of LHS 2397aAB, we 
have determined its total mass to a precision of 10% 
(dominated by the error in the parallax). One way to 
estimate individual masses is to use evolutionary mod- 
els in the same way they were used to infer the sys- 
tem age. By constraining the model-derived individual 
masses of LHS 2397aA and LHS 2397aB to add up to 
the observed total mass, while still matching their ob- 
served luminosities, the Tucson models give masses of 
0.0839lj];gj]?^ Mq and 0.061 Mq for LHS 2397aA 
and LHS 2397aB, respectively. The Lyon models give 
masses of 0.0848;^;|]^i^ M© and 0.060t|];^?^ Mq. We can 
also estimate the mass ratio (q = M-^/Mp^) in this fash- 
ion, and the Tucson models give q = 0.73lo;i2> while the 
Lyon models give q = 0.71'^q'i^. This is among the low- 
est field ultracool binary mass ratios (e.g., sec Burgasser 
et al. 2007) and is consistent with the original estimate 
of 0.76 by Freed et al. (2003) based on their photometry 
and estimated age range (2-12 Gyr). 

We can also estimate the individual masses of 
LHS 2397aAB by subtracting the mass of the stellar 



TABLE 5 

Evolutionary Model-derived Properties 



Property Median 68.3% c.l. 95.4% c.l. 



Tucson Models (Burrows et al. 1997) 

System 



Age (Gyrj-^ 


1.5 


-0.6,4.1 


-0.9,8.5 


q (Mb/Ma) 


0.73 


-0.12,0.16 


-0.19,0.18 


ATeft (K) 


1130 


-40, 30 


-50, 60 


Component A 


Ma (Mq) 


0.0839 


-0.0015, 0.0007 


-0.005, 0.001 


Teff.A (K) 


2580 


-30, 30 


-50,50 


log(ffA) (Cgs) 


5.381 


-0.014, 0.009 


-0.04,0.02 


Ra (Rq) 


0.0978 


-0.0012,0.0012 


-0.002, 0.003 


LiA/Lio 


0.00 


-0.00,0.00 


-0.00, 0.00 


(Jomponent B 


Mb (M0) 


0.061 


-0.011,0.014 


-0.019,0.015 


Teff,B (K) 


1450 


-40, 40 


-80, 70 


log(ffB) (cgs) 


5.35 


-0.13,0.14 


-0.25,0.15 


Rb (Rq) 


0.0868 


-0.005, 0.005 


-0.005, 0.009 


LiB/Lio 


0.5 


-0.5,0.5 


-0.5,0.5 



Lyon Models (DUSTY; Chabrier et al. 2000) 

System 



Age (Gyr)" 


1.8 


-0.8,8.2 


-1.1,8.2 


q (Mb/Ma) 


0.71 


-0.14,0.09 


-0.21,0.10 


AT,,! (K) 


1040 


-40.50 


-50. 80 


txjmpoiioiit .V 


Ma (Mq) 


0.0848 


-0.0012,0.0010 


-0.004, 0.002 


Teff,A (K) 


2470 


-30,30 


-50, 60 


log(ffA) (cgs) 


5.307 


-0.008,0.007 


-0.03,0.01 


Ra (R©) 


0.1073 


-0.0010,0.0011 


-0.002, 0.002 


LiA/Lio 


0.00 


-0.00, 0.00 


-0.00, 0.00 


Component B 


Mb (Mq) 


0.060 


-0.012,0.008 


-0.020, 0.009 


Teff.B (K) 


1430 


-50,40 


-90, 70 


log (3b) (cgs) 


5.31 


-0.15,0.09 


-0.28,0.09 


Rb (R©) 


0.0895 


-0.0041,0.0064 


-0.004,0.012 


Lie /Lio 


0.0 


-0.0,0.9 


-0.0,1.0 



''Both sets of evolutionary models are only computed up 
to an age of 10 Gyr; therefore, this defines the upper limit 
on the model-derived ages. 



primary from the total mass. There are only three 
stars of comparable spectral type and absolute magni- 
tude to LHS 2397aA with dynamical mass estimates: 
LHS 1070B (M8.5, My = 19.24 ± 0.07 mag, Mk = 
10.42 ±0.04 mag; Leinert et al. 1994, 2000), LHS 1070C 
(M9, My = 19.62 ± 0.08 mag, Mk = 10.76 ± 0.04 mag; 
Leinert et al. 1994, 2000), and GJ 1245C (My = 
18.55 ± 0.05 mag, Mk = 9.99 ± 0.04 mag; Henry & 
McCarthy 1993). The spectroscopically unclassified ob- 
ject GJ 1245C is the closest in absolute magnitude to 
LHS 2397aA and the only one of these with a direct 
dynamical mass measurement (0.074±0.013 Mq; Henry 
et al. 1999). LHS 1070BC only has a total mass deter- 
mined from its relative orbit (0.157±0.009 Mq; Seifahrt 
et al. 2008). If we simply assume that the mass of 
LHS 2397aA is the same as GJ 1245C, then the result- 
ing mass estimate for LHS 2397aB is 0.072±0.019 Mq. 
This is somewhat larger than (+0.6a) but consis- 
tent with the mass estimates from evolutionary models 
(O.OeitHn M© and 0.060tHi2 M©), which have some- 
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what smaller uncertainties. This approach is limited by 
the small number and low precision of dynamical masses 
for stars at the bottom of the main sequence, which will 
be improved by future mass measurements for late-M 
dwarf binaries. 

The mass of LHS 2397aB is estimated to be substellar 
as it is below 0.070 0.092 M0, the range in the hydrogen- 
fusing minimum mass for plausible values of the metal- 
licity and helium fraction (Burrows et al. 2001). This is 
not surprising since theoretical models predict that stars 
at the very bottom of the main sequence could not be 
as faint as LHS 2397aB. We can therefore estimate an 
upper limit on the total mass from the hydrogen-fusing 
minimum mass (i=a0.075 Mq for solar metallicity) and 
the highest mass expected for LHS 2397aA. There are 
few mass measurements at the bottom of the main se- 
quence, but we estimate that a conservative upper limit 
on the mass of LHS 2397aA is 0.095 Mq. G1 866C has a 
mass of 0.0930±0.0008 Mq and at My = 17.4±0.4 mag 
is more than a magnitude brighter than LHS 2397aA; no 
objects with masses > 0.100 M© come within 2 mag of 
LHS 2397aA's My. These upper limits on the masses 
of the individual components of LHS 2397aAB imply 
Mtot < 0.170 Mq, 1.7a larger than our derived dynami- 
cal mass. Thus, we suggest that the high-mass tail of the 
Mtot distribution, while formally allowed by the data, is 
practically excluded based on the properties of the indi- 
vidual components. 

4.3. Temperatures and Surface Gravities 

Without radius measurements for LHS 2397aA and 
LHS 2397aB, we cannot directly determine their effec- 
tive temperatures or surface gravities. However, we 
have used evolutionary models to estimate these prop- 
erties in the same fashion as our model-derived age and 
individual masses. 

The Lyon models give effective temperatures for 
LHS 2397aA and LHS 2397aB of 2470±30 K and 
1430±40 K, while the Tucson models give systemat- 
ically hotter but formally consistent temperatures of 
2580±30 K and 1450±40 K. Because of the nearly flat 
mass-radius relation for such low mass objects, it is es- 
sentially the uncertainties in the component luminosi- 
ties that determine the precision in the model-derived 
effective temperatures for LHS 2397aA. Since brown 
dwarfs cool over time, the precision in LHS 2397aB's ef- 
fective temperature also depends somewhat on the un- 
certainty in the age, which in turn depends on the mass 
error. While our uncertainties were derived specifically 
for this system, approximate scaling relations derived by 
Liu et al. (2008) can provide simple estimates of the rel- 

^® If LHS 2397a were a higher order multiple, then its total mass 
would be more consistent with the high-mass tail of the iWtot distri- 
bution. However, the fact that the luminosities and colors of both 
components agree well with single field dwarfs makes unresolved 
multiplicity unlikely. 

Since LHS 2397aA is a chromospherically active star in a bi- 
nary system, it may be possible to estimate its radius using the 
technique employed by Berger et al. (2009) who measured the ro- 
tation period of 2MASSW J0746425+200032A from its chromo- 
sphcric radio emission and combined this with its vsm(i) and or- 
bital inclination (i) to derive its radius. This method assumes 
that the orbital and rotation axes are aligned. With a tisin(j) of 
20 km s-^ (Mohanty & Basri 2003), LHS 2397aA is expected to 
have a rotation period of about 240 minutes. 



ative dependencies of the mass and luminosity errors on 
the uncertainty in Toff. 

The Lyon models give surface gravities for 
LHS 2397aA and LHS 2397aB of log{g) = 5.307^^;^^^ 
and 5.31^Q 51 (cgs), while the Tucson models give sys- 
tematically higher gravities of \og{g) — 5.3811;q q°4 and 

5.35tQ'i[3 (cgs). Since the radius is essentially constant 
with age, the precision in model-inferred surface gravity 
is driven by the precision in the measured total mass. 
The diff'erence between the two sets of model-inferred 
surface gravities, which are formally inconsistent for 
LHS 2397aA, arises from small differences (< 9%) in the 
model-predicted radii (Table 5). 

4.3.1. Comparison to Field Dwarfs 

The effective temperatures we derive from evolutionary 
models for LHS 2397aA and LHS 2397aB can be com- 
pared to those which have been determined for other 
objects of similar spectral type. Such estimates in the 
literature utilize the nearly flat mass radius relationship 
predicted by theoretical models for brown dwarfs, adopt- 
ing either a typical age (e.g., 3 Gyr; Golimowski et al. 
2004) or radius (e.g., 0.90±0.15 Rq; Vrba et al. 2004). 

For LHS 2397aA (M8±0.5), held M7 M9 dwarfs with 
Lboi measurements have estimated effective tempera- 
tures ranging between 1850 2650 K (Leggett et al. 
2001). For LHS 2397aB (L7±l), field L6-L8 dwarfs 
with Lboi measurements have estimated effective temper- 
atures ranging between 1300 1700 K. Both of these broad 
ranges of effective temperature are consistent with our 
mo del- inferred effective temperatures. However, since 
both estimates are based on evolutionary models, this 
only means that field L6-L8 dwarfs from previous studies 
encompass objects of the same mass/age as LHS 2397aB. 

Gautier et al. (2007) used their 24 /xm Spitzer/MIPS 
photometry to determine effective temperatures for 
nine M7.5-M8.5 dwarfs (including LHS 2397a) using 
the nearly model independent infrared flux method. 
These temperatures range from 2400 2730 K (exclud- 
ing LHS 2397a). They derived an effective temperature 
of 2380 K for LHS 2397a, but LHS 2397aB likely affects 
this estimate since it is not a negligible source of flux 
at 24 /xm: the atmospheric models of Burrows et al. 

(2006) predict that it would contribute 0.16 mJy of the 
measured 0.99 mJy at 24 /im. This is likely why the 
Ks - [24] color of LHS 2397a (1.09±0.10 mag) is redder 
than any other object of its spectral type (mean and rms 
of 0.80±0.08 mag), and this could explain why the tem- 
perature estimate for LHS 2397a is outside the range of 
the other eight M7.5 M8.5 dwarfs in the Gautier et al. 

(2007) sample. The infrared flux method temperature 
range of the other M7.5-M8.5 dwarfs (2400-2730 K) is 
in excellent agreement with the evolutionary model in- 
ferred effective temperatures for LHS 2397aA. 

4.3.2. Comparison to Atmospheric Models 

Effective temperatures have also been determined for 
many held dwarfs through spectral synthesis modeling. 
By comparing these results to our evolutionary model- 
derived temperatures, we perform a consistency check 
between these two classes of models, which are essentially 
independent. Even though atmospheric models provide 
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an important boundary condition for evolutionary mod- 
els, very similar bulk properties (e.g., Lboi and Tgs) are 
predicted by evolutionary models using widely varying 
boundary conditions (Chabrier & Baraffe 2000; Saumon 
& Marley 2008). 

Using spectral synthesis modeling over a very nar- 
row spectral range (2.297-2.310 /zm), Jones et al. (2005) 
found effective temperatures of 2900 K for the two M7~ 
M9 dwarfs in their study. Using a broader spectral range 
(1.0-2.5 /im), Leggett et al. (2001) found much cooler 
temperatures of 2100-2300 K for the five M7-M9 dwarfs 
in their study. The evolutionary model-inferred range of 
effective temperatures for LHS 2397aA (2440-2610 K) 
is inconsistent with both results from spectral synthesis 
modeling. Both authors cite the treatment of dust as 
the likely cause of discrepancies they found between the 
atmospheric models and their spectra. 

Using spectral synthesis modeling over a very broad 
spectral range (0.95-14.5 /itm), Gushing et al. (2008) de- 
rived effective temperatures for several field L and T 
dwarfs, including two late-L dwarfs of comparable spec- 
tral type and JHK color to LHS 2397aB. These two ob- 
jects (2MASS J0825-H2115, L7.5; DENIS J0255-4700, 
L8) were both found to have effective temperatures of 
1400 K, which is in perfect agreement with our model- 
derived temperature estimates for LHS 2397aB. This 
is somewhat surprising as previous studies of L and 
T dwarfs have found large (100-200 K) discrepancies be- 
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Fig. 11. — The Hertzsprung-Russell diagram show- 
ing isomass lines from evolutionary models for the mass 
of LHS 2397aB with dotted lines encompassing the la 
uncertainty in the mass. The effective temperature of 
field L7±l dwarfs as determined from spectral synthesis 
(1400 K) is shown as a filled circle with 50 K error bars, 
which correspond to the atmospheric model grid size. 
The open squares demarcate the median and Icr confi- 
dence limits on the evolutionary model-derived age of 
LHS 2397aAB, using the combined constraint of the to- 
tal mass and individual luminosities (see Figure 10). The 
atmospheric model temperature is in very good agree- 
ment with the evolutionary models, unlike previous stud- 
ies which have found 100-200 K discrepancies on the H-R 
diagram (Liu et al. 2008; Dupuy et al. 2009). 



tween temperatures derived from atmospheric and evolu- 
tionary models (Liu et al. 2008; Dupuy et al. 2009). Fig- 
ure 11 shows the atmospheric model predicted tempera- 
ture for LHS 2397aB on the Hertzsprung-Russell (H-R) 
diagram with the evolutionary model tracks. It is inter- 
esting to note that both components of the L4-)-L4 binary 
HD 130948BC were found to lie above the evolutionary 
model tracks (Dupuy et al. 2009), both components of 
the T5-HT5.5 binary 2MASS J1534-2952AB were found 
to lie below the tracks (Liu et al. 2008), and LHS 2397aB 
with an intermediate spectral type of L7 falls exactly on 
the tracks. This suggests that whatever different sys- 
tematic errors are harbored in each class of models may 
cancel out on the H-R diagram at this spectral type. 

Gushing et al. (2008) also derived surface gravities for 
the two late-L dwarfs in their study by direct model fit- 
ting {\og{g) — 4.5-5.5) and by using their evolutionary 
sequences (\og{g) = 5.4-5.5). The latter values were de- 
termined from the atmospheric model T^g and the radius 
derived from the flux normalization constant {R/d)"^, 
which yields the radius if the distance is known. This 
is because atmospheric models predict the emergent flux 
in absolute units, which can be compared directly to flux- 
calibrated spectra. With independent estimates of both 
R and Tcff (from atmosphere fitting), any other property 
can be derived using evolutionary models. The surface 
gravities derived by Gushing et al. (2008) in this way for 
these two objects are higher by 0.4-1.9cr than our evolu- 
tionary model inferred values for LHS 2397aB from Lyon 
and Tucson models {\og{g) = 5.31-5.34). 

4.4. Near- Infrared Colors 

The Lyon evolutionary models provide predictions of 
the fluxes of ultracool dwarfs in standard filter band- 
passes as a function of model mass and age.^^ We have 
derived the model-predicted near-infrared colors of both 
components of LHS 2397aAB in the same fashion as the 
individual masses, effective temperatures, and age (i.e., 
using the combined observational constraints of the to- 
tal mass and individual luminosities). For LHS 2397aA, 
the Lyon DUSTY models predict a very narrow range of 
colors: J - K = 0.784 ± 0.007 mag, H - K = 0.286 ± 
0.002 mag, and J-H = 0.498±0.006 mag. None of these 
are consistent with the observed colors of LHS 2397aA: 
J - K ^ 1.28 ± 0.03 mag, H - K = 0.43 ± 0.03 mag, 
and J — H = 0.85 ± 0.03 mag. All of the observed colors 
are redder than predicted, with J — K being the most 
discrepant (0.50 mag), followed by J — H (0.35 mag) and 
H-K (0.14 mag). Since the J-fsT color of LHS 2397aA 
is consistent with other field dwarfs of spectral type M8 
(1.14±0.10 mag; e.g., Faherty et al. 2009), the Lyon 
DUSTY models will generally provide inaccurate esti- 
mates of the fundamental properties of late-M dwarfs 
from their near-infrared colors. In fact, it would not be 
possible to estimate the mass of LHS 2397aA from the 
color-magnitude diagram as no model isochrone is con- 
sistent with the data. 

We show the model-predicted JHKL' colors of 

The models give JHK photometry on the CIT system, and 
we convert our measured MKO system colors to this system using 
the relations of (Stephens & Leggett 2004). The model L'-band 
photometry is given on the Johnson-Glass system, and no empirical 
relations between this and the MKO photometric system exist, so 
we neglect any differences. 
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Fig. 12. — Color-magnitude diagrams showing the measm'ed properties of LHS 2397aB compared to Lyon evolutionary 
tracks. The solid lines are isomass tracks from the DUSTY (Chabrier et al. 2000) and COND (Baraffe et al. 2003) 
models [JHK photometry on the CIT system and i'-band photometry on the Johnson-Glass system). The open 
squares demarcate ages along the isomass lines, showing the median and la confidence limits on the model-inferred 
age of LHS 2397aAB if possible (0.5 Gyr also shown for reference). Field dwarfs with parallax measurements better 
than 10% and JHK colors better than 0.10 mag are shown as filled gray circles. LHS 2397aB has near-infrared colors 
typical of field dwarfs, but neither set of evolutionary models reproduces the observed colors of LHS 2397aB. This is 
not surprising since each model is a limiting case in the amount of dust in the photosphere (very dusty or no dust) and 
LHS 2397aB is expected to be intermediate given its estimated spectral type of L7±l. The K — L' colors predicted 
by models arc by far the least discrepant with the data. 



LHS 2397aB compared to the observations and other 
field dwarfs in Figure 12. The colors of LHS 2397aB 
are fully consistent with field dwarfs of similar i^-band 
absolute magnitude but are significantly bluer than the 
DUSTY models for any assumed age. This discrepancy 
with the models is not surprising as it is well known 
that the DUSTY models do not accurately reproduce 
the near-infrared colors of late-L dwarfs. This is at least 
partly because DUSTY models represent the limiting 
case of very high dust content in the photosphere, while 
late-L dwarfs are actually in transition to dust-free at- 
mospheres. For comparison, we also show the COND 
models (Baraffe et al. 2003), which represent the con- 
verse limiting case of a dust-free photosphere. As such, 
the COND models reproduce the colors of the bluest 
T dwarfs well, in which dust has sedimented below the 
photosphere, but they fail to reproduce the colors of 
LHS 2397aB. For both DUSTY and COND models, the 
least discrepant color is K — L' , which is only about 
0.5 mag bluer than predicted by DUSTY or COND, pre- 
sumably because the effects of dust are less severe at 
longer wavelengths. The difficulty in modeling the colors 
of objects like LHS 2397aB that are in transition from 
the L to T spectral class is well known, though mod- 
els developed recently that attempt to account for the 
sedimentation of dust have enjoyed some success at re- 
producing observations (Burrows et al. 2006; Saumon & 
Marley 2008). 



5. DISCUSSION 



5.1. 



Effective Temperature of the L/T Transition 

Measurements of the fundamental properties of brown 
dwarfs undergoing the transition from L to T dwarfs are 
critical for understanding the physics governing dust in 
ultracool atmospheres. With a spectral type of L7±l, 
LHS 2397aB is one such object for which our mass and 
luminosity measurements have enabled precise effective 
temperature and surface gravity estimates (Section 4.3). 
In fact, Liu et al. (2008) noted that "mass benchmarks" 
such as LHS 2397aB provide much more precise tem- 



perature estimates than "age benchmarks" (i.e., substel- 
lar companions to stars of known age) since the ages of 
field stars are notoriously difficult to estimate precisely 
(e.g., see Mamajek & Hillenbrand 2008, and references 
therein). Metchev & Hillenbrand (2006) suggested that 
the effective temperature of the L/T transition varies 
with surface gravity (or equivalently, age) , based on such 
age benchmark objects. Thus, it is interesting to com- 
pare the estimated temperature of LHS 2397aB to those 
of the few late-L and early-T dwarf companions with 
age estimates: Gl 584C (L8; Kirkpatrick et al. 2001), 
HD 203030B (L7.5; Metchev & Hillenbrand 2006), and 
HN Peg B (T2.5; Luhman et al. 2007). 

For Gl 584C, Golimowski et al. (2004) derived 
an effective temperature range of 1300-1400 K from 
evolutionary models, using its measured luminosity 
(logLboi/Lo = — 4.58±0.04) and estimated age range (1- 
2.5 Gyr). This is consistent with our temperature esti- 
mate for LHS 2397aB (1380-1470 K from Lyon models; 
1410-1490 K from Tucson models). In a similar fashion, 
Metchev & Hillenbrand (2006) and Luhman et al. (2007) 
estimated the effective temperature ranges for two com- 
panions to young stars: 1090-1280 K for HD 203030B 
(130-400 Myr) and 1060-1200 K for HN Peg B (100- 
500 Myr). Compared to Gl 584C and LHS 2397aB, 
the temperatures of these young companions are much 
cooler. As discussed by Metchev & Hillenbrand (2006), 
field dwarfs of similar spectral type (L7-L9) with mea- 
sured luminosities also have higher estimated tempera- 
tures (mean and rms of 1460±100 K; from Vrba et al. 
2004), though these estimates depend strongly on the 
estimated age of the field population, which is not well- 
constrained observationally. 

Metchev & Hillenbrand (2006) suggested two possi- 
ble explanations for the observed differences in effective 
temperatures: (1) the age of the field population has 

19 Gl 337CD (L8+T; Wilson et al. 2001; Burgasser et al. 2005) 
is another such companion; however, it is a tight binary with poor 
constraints on its resolved properties. Thus, we exclude it from our 
discussion. 
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been overestimated and (2) the effective temperature of 
the L/T transition depends on surface gravity, ff field 
dwarfs are younger than estimated, their radii will be 
larger and thus lower effective temperatures are required 
to match their measured luminosities. Support for this 
possibility has come from Liu et al. (2008), who sug- 
gested that the age of field dwarfs may be overestimated 
by as much as a factor of 6±3 (implying an age range 
of 0.3-1.0 Gyr for field objects) based on a wlOO K dis- 
crepancy between the field population and their dynam- 
ical mass-based temperature estimates for the compo- 
nents of 2MASS J1534-2952AB (T5+T5.5). Such an 
overestimate would impact effective temperature esti- 
mates at all spectral types, not just the L/T transition. 
Metchev & Hillenbrand (2006) preferred the explanation 
that the L/T transition is gravity dependent, so that 
objects of different radii/ages would undergo this evolu- 
tionary phase at different effective temperatures. Our ef- 
fective temperature determination for LHS 2397aB sup- 
ports this idea since LHS 2397aB belongs to the field 
population (Section 3.4) and its temperature is incon- 
sistent by 100 400 K with the young late-L companions 
HD 2030306 and HN Peg B. 

5.2. Direct Measurement of the Mass Ratio 

Although the total dynamical mass we have measured 
enables strong tests of theoretical models, these tests 
would be even more potent with the direct measurement 
of the individual component masses. In principle, this 
can be done using radial velocity or astrometric moni- 
toring, both of which provide a measurement of the bi- 
nary's mass ratio. The radial velocity semi-amplitude of 
LHS 2397aA is 1.9 km s^^ for the mass ratio of 0.7 de- 
rived from evolutionary models (Section 4.2). Thus, a 
radial velocity precision of 0.1 km s~^ would be needed 
to achieve a mass ratio measurement precise to 10% (i.e., 
comparable to the precision in the total mass). As dis- 
cussed in Section 3.1, Basri & Reiners (2006) measured 
radial velocities for LHS 2397aA at two epochs separated 
by almost exactly half an orbital period, but with a pre- 
cision of 1.6 km s~^ these data are not sufficient to de- 
termine the mass ratio. 

To determine the individual masses from astrometric 
monitoring, the orbit of LHS 2397aA about the system's 
barycenter must be measured. The size of this orbit is 
89 mas for a mass ratio of 0.7. Therefore, an astrometric 
precision of 5 mas in the optical, where LHS 2397aB is es- 
sentially invisible, is all that would be needed to measure 
the mass ratio to «10% {(^q/q = (1 + q)(Jai/cii)- Such 
a measurement is quite feasible, though it requires long- 
term observations given the orbital period of 14.2 years. 

5.3. Lithium Depletion 

The presence or absence of lithium has long been used 
as a tracer of the internal structure of objects, provid- 
ing constraints on mass and/or age, because it is only 
destroyed by fusion at temperatures above f«2.5xl0® K. 
There is a mass limit below which lithium is never de- 
stroyed in the fully convective interiors of very low-mass 
stars, and this occurs at «0.06 Mq (Chabrier et al. 1996). 
The estimated mass of LHS 2397aB is extremely close to 
this limit, so Icr confidence limits on its model-derived 
lithium abundance given in Table 5 range from zero (fully 
depleted) to unity (no depletion). There are currently 



no objects with direct mass estimates to provide obser- 
vational constraints on the lithium depletion boundary 
for brown dwarfs. Thus, observations of lithium absorp- 
tion at 6708 A for LHS 2397aB would provide a strong 
test of theoretical models of brown dwarf interiors. In 
addition, such observations would provide an indepen- 
dent constraint on the mass and age of LHS 2397aB that 
would enable a test of model-predicted luminosity evo- 
lution. Lithium observations are only possible with a 
repaired Space Telescope Imaging Spectrograph (STIS) 
onboard HST, and the next few years are an ideal time 
for such a measurement as the orbital separation is cur- 
rently approaching its maximum (Figure 6). 

6. CONCLUSIONS 

We have determined the orbit of the M8-I-L7 bi- 
nary LHS 2397aAB using relative astrometry spanning 
11.8 years of its 14.2-year orbit. The astrometry and 
corresponding errors used to derive this orbit were thor- 
oughly examined through Monte Carlo simulations. The 
resulting best-fit orbit has a reduced of 1-04 and to- 
tal mass of 0.146^0 013 The error in the dynamical 
mass (10%) is dominated by the 3.0% error in the paral- 
lax, which translates into a 9.0% error in the mass. The 
total mass is consistent with the primary component be- 
ing a star and the companion being a brown dwarf, as 
expected from their spectral types. The combined obser- 
vational constraints of the total mass and individual lu- 
minosities break the mass age - luminosity degeneracy for 
the substellar companion, enabling specific predictions of 
its properties from evolutionary models. Tucson models 
(Burrows et al. 1997) predict an age for the system of 
l-5lo:6 Gyr, and the Lyon (DUSTY; Chabrier et al. 2000) 
model-derived age of 1.81q;| Gyr is consistent. These 
ages are consistent with the (weak) observational con- 
straints on the system: LHS 2397a's space motion im- 
plies it is a member of the thin disk (1-10 Gyr; e.g., 
Wood & Oswalt 1998), and although it is chromospher- 
ically active, its expected activity lifetime is > 8 Gyr 
(West et al. 2008). 

LHS 2397aB (L7±l) is now the first mass benchmark 
at the transition from L to T spectral types. This enables 
a precise estimate of the effective temperature of this 
phase of substellar evolution from theoretical models, 
which give 1450±40 K (Tucson) and 1430±40 K (Lyon). 
This is w200 K higher than temperature estimates for 
late-L companions to young stars: HD 203030B (130- 
400 Myr; Metchev & Hillenbrand 2006) and HN Peg B 
(100-500 Myr; Luhman et al. 2007). However, the ef- 
fective temperature of LHS 2397aB is consistent with 
the 1300-1400 K range estimated for the older late-L 
companion Gl 584C (1-2.5 Gyr; Kirkpatrick et al. 2001). 
This supports the idea originally proposed by Metchev 
& Hillenbrand (2006) that the L/T transition can oc- 
cur over a wide range of effective temperatures due to a 
surface gravity dependence on the process. In addition, 
the temperature for LHS 2397aB is consistent with esti- 
mates for field late-L dwarfs (1440±100 K), and thus an 
overestimate of the age of the field population (or under- 
estimate of radii) is not required to explain the range of 
temperatures observed at the L/T transition, as previ- 
ously suggested (Metchev & Hillenbrand 2006; Liu et al. 
2008). 
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Wc have determined the spectral type of LHS 2397aB 
from its JHKL' absolute magnitudes, using a novel tech- 
nique that enables a quantitative assessment of the spec- 
tral type uncertainty by accounting for both measure- 
ment errors and the rms scatter in the empirical rela- 
tions. Comparing to other objects of similar spectral 
type and color for which Gushing et al. (2008) have con- 
ducted spectral synthesis analysis allows us to test atmo- 
spheric models for consistency with evolutionary mod- 
els. The Gushing et al. (2008) study included two late- 
L dwarfs, and their derived effective temperatures were 
both 1400 K, which is consistent with the effective tem- 
perature of LHS 2397aB we derived. This is somewhat 
surprising as similar consistency tests by previous studies 
of binary brown dwarfs have found large (100-200 K) dis- 
crepancies between effective temperatures derived from 
atmospheric and evolutionary models. However, it is in- 
teresting to note that LHS 2397aB has a spectral type in- 
termediate between HD 130948BG for which the effective 
temperatures are underestimated by atmospheric models 
(L4-I-L4; Dupuy et al. 2009) and 2MASS J1534-2952AB 
for which the temperatures are overestimated (T5+T5.5; 
Liu et al. 2008). This suggests that LHS 2397aB may be 
at a temperature where atmospheric and evolutionary 
models are in agreement, because the significant system- 
atic errors in the models cancel out. 

Because of its extreme constrast ratio (> 100:1 in the 
optical), LHS 2397aAB is one of the few ultracool bina- 
ries whose mass ratio can be readily measured through 
astrometric monitoring. Given the estimated mass ratio 
of 0.7, the barycentric orbit of LHS 2397aA is 89 mas, 
which is easily mcasureable with current instrumenta- 
tion. It is unlikely that LHS 2397aA's 14-year orbital 
motion significantly impacted the original optical paral- 
lax measurement, except by introducing additional scat- 
ter in the measurements. Thus, new astrometric mon- 
itoring observations should improve the parallax preci- 
sion, which is essential for further testing of models. 
The total mass is currently constrained to a precision 
of 10%, dominated by the error in the parallax, and 
at this precision the secondary component is formally 
allowed (at l-2a) to be a star at the bottom of the 
main sequence. A more precise mass measurement for 
this system would enable miich stronger constraints on 
all model-predicted properties, including the amount of 
lithium depletion in LHS 2397aB. Its model-derived in- 
dividual mass («0.06±0.01 Mq) is very close to the the- 
oretical mass-limit for lithium burning, which has never 
been tested with direct mass measurements. A lithium 
measurement for a brown dwarf of known mass would en- 
able a unique constraint on substellar models as lithium 
depletion proceeds rapidly and independently of lumi- 
nosity evolution. Such a resolved measurement of the 
the lithium doublet at 6708 A for LHS 2397aB is only 
possible with HST/STIS. 
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